Breakdown of the Migdal-Eliashberg theory in the strong-coupling adiabatic regime 
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In view of some recent works on the role of vertex cor- 
rections in the electron-phonon system we readress an impor- 
tant question of the validity of the Migdal-Eliashberg theory. 
Based on the solution of the Holstein model and 1/A strong- 
coupling expansion, we argue that the standard Feynman- 
Dyson perturbation theory by Migdal and Eliashberg with or 
without vertex corrections cannot be applied if the electron- 
phonon coupling is strong (A > 1) at any ratio of the phonon, 
(j and Fermi, Ef energies. In the extreme adiabatic limit 
(w << Ef) of the Holstein model electrons collapse into 
self-trapped small polarons or bipolarons due to spontaneous 
translational-symmetry breaking at A ~ 1. With the increas- 
ing phonon frequency the region of the applicability of the 
theory shrinks to lower values of A < 1. 

PACS numbers:74.20.-z, 74.65. +n,74.60.Mj 



The theory of ordinary metals is based on 'Migdal's' the- 
orem [jjj , which showed that the contribution of the dia- 
grams with 'crossing' phonon lines (so called 'vertex' cor- 
rections) is small. This is true if the parameter Xu>/Ep 
is small. Neglecting the vertex corrections Migdal calcu- 
lated the renormalised electron mass as m* = m(l + A) 
(near the Fermi level) jlj] and by breaking the gauge sym- 
metry Eliashberg extended the Migdal theory to describe 
the BCS superconducting state at intermediate values of 
A 0. The same theory, applied to phonons, yields the 
renormalised phonon frequency Co = uj(1 — 2A) 1 / 2 Jl|] with 
an instability at A = 0.5. Because of this instability both 
Migdal 0] and Eliashberg [|| restricted the applicability 
of their approach to A < 1. 

It was then shown that if the adiabatic Born- 
Oppenheimer approach is properly applied to a metal, 
there is only negligible renormalisation of the phonon fre- 
quencies of the order of the adiabatic ratio, LojEp << 1 
for any value of A ||. The conclusion was that the 
standard electron-phonon interaction could be applied to 
electrons for any value of A but it should not be applied 
to renormalise phonons Q]. As a result, many authors 
used the Migdal-Eliashberg theory with A much larger 
than 1 (see, for example, Ref. ||). 

However, starting from the infinite coupling limit, A = 
00 and applying the inverse (1/A) expansion technique || 
we showed 0-^] that the many-electron system collapses 
into small polaron regime at A ~ 1 almost independent 
of the adiabatic ratio. This regime is beyond the Migdal- 
Eliashberg theory. It cannot be reached by summation 
of the standard Feynman-Dyson perturbation diagrams 
even including all vertex corrections, because of the bro- 
ken translational symmetry, as first discussed by Landau 



fllQf for a single electron and by us |y for the many- 
electron system. During last years quite a few numeri- 
cal and analytical studies have confirmed this conclusion 
pl| pq| (and references theirein). On the other hand a 
few others (see, for example, |29|]3C| 1) still argue that the 
breakdown of the Migdal-Eliashberg theory might hap- 
pen only at A > Ep/w >> 1. Indeed numerical study of 
the finite bandwidth effects |3lJ and some analytical cal- 
culations of the vertex corrections to the vertex function 
P4B3 30 j with the standard Feynman-Dyson perturba- 



tion technique confirm the second conclusion. 

In this letter I compare the Migdal solution of the Hol- 
stein Hamiltonian with the exact one in the extreme adi- 
abatic regime, u/Ep — > 0, to show that the ground state 
of the system is a self-trapped insulating state with bro- 
ken translational symmetry already at A > 1. 

The vertex corrections and the finite bandwidth are 
rather technical issues, playing no role in the extreme 
adiabatic limit 0|32|jl|. There is, however, another ba- 
sic assumption of the canonical Migdal-Eliashberg the- 
ory. That is the electron and phonon Green functions 
(GF) are translationally invariant. As a result one takes 
G(r, r', t) = G(r — r', r) with Fourier component G(k, Q) 
prior to solving the Dyson equations. This assumption 
excludes the possibility of local violation of the transla- 
tional symmetry due to the lattice deformation in any 
order of the Feynman-Dyson perturbation theory. This 
is similar to the absence of the anomalous (Bogoliubov) 
averages in any order of the perturbation theory. To en- 
able electrons to relax into the lowest polaronic states, 
one has to introduce an infinitesimal translationally non- 
invariant potential, which should be set equal to zero 
only in the final solution for the GF JsJ . As in the case 
of the off-diagonal superconducting order parameter, a 
small translational-symmetry-breaking potential drives 
the system into a new ground state at sufficiently large 
coupling, A ~ 1, independent of the adiabatic ratio. Set- 
ting it equal to zero in the solution of the equations of 
motion restores the translational symmetry but in a new 
polaronic band rather than in the bare electron band, 
which turns out to be an excited state. 

In particular, let us consider the extreme adiabatic 
limit of the Holstein chain |34|| , which has the simple 
analytical solution: 



H =-tJ2 4 c i + H - c - + 2(Afci) 1/2 

<ij> i 
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where t is the nearest neighbour hopping integral, c|,Cj 
are the electron operators, Xi is the normal coordinate of 
the molecule (site) i, and k = Mui 2 with M the ion mass 
p5| . We first consider the two-site case (zero dimensional 
limit), i,j = 1,2 with one electron, and than generalise 
the result for the infinite lattice with many electrons. 
The transformation X — {x\ + X2), £ = x\ — x 2 allows 
us to eliminate the coordinate X, which is coupled only 
with the total density (m + n 2 = 1), leaving the following 
Hamiltonian to solve in the extreme adiabatic limit M — > 
00: 



Ck 



H = -t{c\c 2 + eld) + {\kt) l ' 2 ^{c\c x - c\c 2 ) + ^- 



The solution is 



^ = (ac t 1 +/3 C t)|0>, 



where 



[t 2 + ((\kty/ 2 z + (t 2 + xktey/ 2 ) 2 ] 1 ' 2 ' 

{xkt) 1 ' 2 ^ + (t 2 + \ktey' 2 
( ~ ~ [t 2 + ((\kty/ 2 £_ + {t 2 + xktey/ 2 ) 2 ] 112 ' 



and the energy 



E = ^-(t 2 + xktei 



2U/2 



(2) 
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In the extreme adiabatic limit the displacement £ is 
classical, so the ground state energy, Eq and the ground 
state displacement £0 are obtained by minimising Eq.(6) 
with respect to £. If A > 0.5 one obtains 



E 



and 



t(4X 2 - 1) 



-,1/2 
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The symmetry-breaking ('order') parameter is 



A = f3 2 



[2A+(4A 2 -l) 1 / 2 ] 2 -l 
[2A+(4A 2 -l) 1 /2]2 + 1 - 



(7) 
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If, however, A < 0.5 the ground state is translationally 
invariant with Eg = —t, £ = 0, (3 = —a, and A = 0. Pre- 
cisely this state is the 'Migdal' solution of the Holstein 
model. Indeed, in the Migdal approximation GF is diag- 
onal in the k representation, G(k, k',r) = G(k, r)i5k.k' • 
The site operators can be transformed into momentum 
space as 



N x l 2 cj exp(ikaj), 



(10) 



with a the lattice constant, N the number of sites, and 
k = 2-rm/Na with -N/2 < n < N/2. Than the off- 
diagonal GF with k = and k' = 11 /a of the two-site 
chain (N = 2) at r = —0 is given by 

G(fc,fc',-0) = ^<( c t 1 -4)( Cl + C2 )). (11) 
Calculating this average one obtains 



G(k,k',-0) = -(a 2 -f3 2 ), 



(12) 



which should vanish in the Migdal theory. Hence, this 
theory only provides symmetric (translationally invari- 
ant) solution with \a\ — \(3\. When A > 0.5 this solution 
is not the ground state of the system, Fig.l. The system 
collapses into a localised adiabatic polaron, trapped on 
the 'right' (or on the 'left') site due to a finite local lattice 
deformation £o- On the other hand, when A < 0.5, the 
Migdal solution is the only solution with £0 = 0. Thus 
the Migdal-Eliashberg constraint [|]J|] on the applicabil- 
ity of their approach is perfectly correct irrespective of 
the phonon frequency renormalisation. 

The generalisation to the multi-polaron system on the 
infinite lattice of any dimension is straightforward in the 
extreme adiabatic regime. The adiabatic solution of the 
infinite one-dimensional (ID) chain with one electron was 
obtained by Rashba |3q| in the continuous approxima- 
tion, and by Holstein [B4j and Kabanov and Mashtakov 
Jll| for a discrete lattice. The last authors also studied 
the Holstein two-dimensional (2D) and three-dimensional 
(3D) lattices in the adiabatic limit. According to Ref. 
Jilt the self-trapping of a single electron occures for any 
value of A in a ID Holstein chain, and at A > 0.875 and 
A > 0.92 in 2D and 3D, respectively. The single-polaron 
GF is not translationally invariant in this strong-coupling 
adiabatic limit. The radius of the self-trapped adiabatic 
polaron, rv, is readily derived from its continous wave 
function [Bql 



ip(x) ~ 1/ cosh(Ax/a). 



(13) 



It becomes smaller than the lattice constant, r p = a/X for 
A > 1. That is why in the strong-coupling (A > 1) adi- 
abatic regime the multi-polaron system remains in the 
self-trapped insulating state no matter how many po- 
larons it has. The only instability which might occur in 
this regime is the formation of on-site self-trapped bipo- 
larons, if the Holstein on-site attractive interaction, 2Xzt, 
is larger than the repulsive Hubbard U |37[]. Actually, 
this instability can be seen in the second order Feynamn- 
Dyson diagramm containing the polarisation loop, as ex- 
plained in Ref. [3^ ]. On-site bipolarons forme charge 
ordered insulating state due to weak repulsion between 
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them []39|. The exact analytical and numerical proof of 
this statement as well as different polaronic and bipola- 
ronic configurations in the adiabatic Holstein model was 
reviewed by Aubry | [40| |. For example, the asymptoti- 
cally exact many-particle ground state of the half-filled 
Holstein model in the strong-coupling limit (A >> f ) is 

*=iRt4j°> ( m ) 

jeB 

for any value of the adiabatic ratio, uj/zt ||{],|l|. Here j 
are -B-sites of the bipartite lattice A + B. It is a charged 
ordered insulating state, rather than the Fermi liquid, ex- 
pected in the Migdal approximation at any value of A in 
the extreme adiabatic limit, uj — > 0. If the Coulomb re- 
pulsion (Hubbard U) is sufficiently strong to prevent the 
bipolaron formation, than at half filling every site is oc- 
cupied by one polaron. The deformation barrier for their 
tunneling to a neighboring site disappears, so one could 
erroneously believe that the system should be metallic. 
However, to keep polarons unbound into on-site bipo- 
larons, the Hubbard U should be larger than 2Xzt, i.e. 
larger than the bandwidth, U > 2zt, if A > 1. It is well 
known that in this regime the system is a Mott-Hubbard 
insulator, rather than an ordinary metal. Hence, the 
Migdal-Eliashberg theory cannot be applied at A > I, no 
matter what a strength of the Coulomb interaction, the 
number of electrons and the dimension of the lattice are. 

The non-adiabtic corrections (phonons) allow polarons 
and bipolarons propagate as Bloch waves in a new (nar- 
row) band. Thus, under certain conditions |2§f| the multi- 
polaron system might be metallic with polaronic (or bipo- 
laronic) carriers rather than bare electrons. However, 
there is a qualitative difference between the ordinary 
Fermi liquid and the polaronic one. In particular, the 
renormalized (effective) mass of electrons is independent 
of the ion mass M in ordinary metals (where the Migdal 
adiabatic approximation is believed to be valid), because 
A does not depend on the isotope mass. However, the 
polaron effective mass m* will depend on M. This is be- 
cause the polaron mass m* — mexp(A/uj) where m 
is the band mass in the absence of the electron-phonon 
interaction, and A is a constant. Hence, there is a large 
isotope effect on the carrier mass in polaronic metals, 
in contrast to the zero isotope effect in ordinary metals. 
Recently, this effect has been experimentally found in 
cuprates Q and manganites Q. With the increasing 
phonon frequency the 1/A polaron expansion becomes 
valid for a smaller values of A (see, for example, [p8[). 
Hence, the region of the applicability of the Migdal ap- 
proach shrinks with increasing u> to the smaller values of 
the coupling, A < I. 

The essential physics of strongly coupled electrons and 
phonons has been understood with the 1/ A expansion 
technique p,^6|, which starts with the exact solution in 
the extreme limit A = oo and allows for the summation 



of all multiphonon diagrams in any order of the small 
parameter I/A for any value of the adiabatic ratio uj/t. 
It predicts a breakdown of the Migdal-Eliashberg the- 
ory at A ~ 1 or less, in agreement with the exact solu- 
tion of the extreme adiabatic Holstein model discussed 
here. There are also other studies of the same problem, 
which do not relay on the standard Feynman-Dyson per- 
turbation theory in powers of A. In particular, Takada 
et al |T^j2^,Q applied the gauge-invariant self-consistent 
method neglecting the vector term in the Ward identity. 
Benedetti and Zeyher ]23| ] applied the dynamical mean- 
field theory in infinite dimensions. As in the I/A expan- 
sion technique, both approaches avoided the problem of 
the broken translational symmetry by using the nondis- 
persive vertex and Green functions as the starting point. 
As a result they arrived at the same correct conclusion 
about the applicability of the Migdal approach (in Ref. 
p3| the critical value of A was found to be 1.3 in the 
adiabatic limit). Meanwhile, some authors |30| , p9|| , who 
relay on the perturbation expansion in powers of A, fail 
to notice the self-trapping transition, erroneously claim- 
ing that the Migdal-Eliashberg theory could be applied 
at practically any A in the extreme adiabatic regime. 

In conclusion, I have shown that the Migdal-Eliashberg 
theory cannot be applied in the adiabatic regime of the 
strongly coupled many-electron system if the (BCS) cou- 
pling constant, A, is about 1 or larger irrespective of the 
number of electrons. There is a transition into the self- 
trapped state due to the broken translational symmetry, 
Fig.l. The transition appears at 0.5 < A < 1.3 depending 
on the lattice dimensionality. Naturally, with the increas- 
ing phonon frequency the region of the applicability of 
the Migdal approach shrinks further to lower values of A. 
I believe that this conclusion is correct for any electron- 
phonon interaction conserving the on-site electron occu- 
pation numbers. In particular, Hiramoto and Toyozawa 
|15| calculated the strength of the deformation poten- 
tial, which transforms electrons into small polarons and 
bipolarons. Their continuous approach is sufficient for a 
qualitative estimate if the Debye wavenumber qjj ~ 7r/a 
is introduced as an upper limit cut-off in all sums in the 
momentum space. Hiramoto and Toyozawa found that 
the transition between two electrons and a small bipo- 
laron occurs at A ~ 0.5, that is half of the critical value of 
A at which the transition from electron to small polaron 
takes place in the extreme adiabatic limit (sqo << zt, 
s the sound velocity). The effect of the adiabatic ratio 
sqo/zt on the critical value of A was found to be negligi- 
ble. The radius of the acoustic polaron and bipolaron is 
the lattice constant, so the critical A docs not depend on 
the number of electrons in this case either. 

The author greately appreciates enlightening discus- 
sions with A.R. Bishop, V.V. Kabanov, R. Zeyher, and 
J.R. Schrieffer. 
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Figure caption 

Fig.l. The ground state energy (solid line), E/t, and 
the 'order' parameter (thin line), A, of the adiabatic 
Holstein model as functions of the coupling constant, A. 
There is the symmetry breaking transition at A = 0.5. 
The energy of the symmetric state is shown by the dashed 
line. 
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